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A cavity perturbation technique is used to study the microwave response 
of the organic superconductor n— (BEDT-TTF) 2 Cu(NCS) 2 . Observation of a 
Josephson plasma resonance, below T c (~10 K), enables investigation of the 
vortex structure within the mixed state of this highly anisotropic, type— II, su- 
perconductor. Contrary to previous assumptions, frequency dependent stud- 
ies (28 — 153 GHz) indicate that the squared plasma frequency {(Ji) depends 
exponentially on the magnetic field strength. Such behavior has been pre- 
dicted for a weakly pinned quasi-two-dimensional vortex lattice [Bulaevskii et 
al. Phys. Rev. Lett. 74, 801 (1995)], but has not so far been observed ex- 
perimentally. Our data also suggests a transition in the vortex structure near 
the irreversibility line not previously reported for an organic superconductor 
using this technique. 

PACS numbers: 71.18.+y, 71.27.+a, 74.25.Nf 

It is widely acknowledged in the layered high— T c superconductors (HTSs) and, more 
recently, in the layered organic superconductors (OSs), that Josephson coupling is respon- 
sible for interlayer transport of superconducting pairs |l],|2].|3|Jl],[|] . Magneto-optical exper- 
iments have been carried out by many groups, and for a wide range of materials, e.g. 
Bi 2 Sr 2 CaCu 2 8+(5 (BSCCO) ||J YBa 2 Cu 3 6 _j / (YBCO) [§, as well as several OSs §§. In 
each case, these experiments have demonstrated that a Josephson plasma resonance ( JPR) is 
observed in the inter-layer conductivity, and that this resonance may be used as an extremely 
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sensitive tool for probing the vortex structure/dynamics. 

/c-(BEDT-TTF) 2 Cu(NCS) 2 (BEDT-TTF denotes bis-ethylenedithio-tetrathiafulvalene, 
or ET for short) is one in a series of highly anisotropic OSs, characterized by a layered 
structure. The conducting layers consist of a checkerboard pattern of face-to-face ET dimers, 
with nearest neighbor dimer pairs oriented perpendicularly |J. These donor sheets exhibit 
near isotropic 2D conductivity within the be plane. The insulating layers form from weakly 
bonded arrays of V— shaped Cu(NCS)2 anions. The anisotropy parameter in the normal 
state, given by the ratio of the in-plane to out-of-plane conductivities <7{, c /cr a , is ~ 1000. In 
the superconducting state, the anisotropy parameter is given by 7 = Aj_/A|| ~ 100 — 200, 
where An and Aj_ are the London penetration depths for AC currents induced parallel (An ~ 
0.8/xm) and perpendicular (Aj_ ~ 100//m) to the conducting layers respectively 0. Such a 
large anisotropy makes this OS a prime candidate to study Josephson coupling, and changes 
in this coupling upon the introduction of supercurrent vortices into the sample through the 
application of an external magnetic field. 

In highly anisotropic superconductors, it is expected that the plasma mode for the low 
conductivity direction (a a ) lies below the superconducting gap i.e. hu p < 2A. Below the 
critical temperature T c , this plasma mode dominates the a-axis microwave response with 
frequency, u p , which depends on the maximum inter-layer (or Josephson) current J m (B,T), 
through the expression 

^(B,T) = ^J m (S,T), 



where ( 1 ) 

Jm = J ((COS (r)) t )rf, 

s is the crystal inter-layer spacing, e c is the high frequency permittivity, $ is the flux 
quantum, <£>„ iri+ i(r) is the gauge-invariant phase difference between layers n and n + 1 at a 
point r = x,y in the be (high conductivity) plane, and (....)* and (....)<* denote thermal and 
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disorder averages. J (T) = c$ /87r 2 sA^(T) is the maximum inter-layer Josephson current 
density at zero field (B DC = 0), and Aj_ is the inter-layer London penetration depth 

An important property of the Josephson coupling is the influence of an externally applied 
magnetic field on the collective plasma oscillation frequency uj p . When the applied DC 
magnetic field is parallel to the least conducting axis, a mixed state is created in which 
the field penetrates the sample in quantized flux tubes, generating supercurrent vortices in 
the superconducting layers. <y9 n n+ i(r) depends explicitly on the vortex structure within this 
mixed state and is, thus, responsible for the field dependence of the resonance frequency uo p 
||. If the flux tubes form straight lines along the a-axis, (cos<£> n . n +i(r)) = 1, and maximum 
Josephson coupling occurs. However, in the presence of disorder, e.g. as a result of crystal 
defects which create vortex pinning sites, or through thermal fluctuations, the flux tubes 
will deviate from straight lines. This suppresses the maximum Josephson current. It is when 
AC currents are excited between the layers, at a frequency which corresponds to the natural 
frequency of the plasma oscillation (u p ), that a sharp resonance is observed. Consequently, 
this resonance frequency provides a direct measurement of the maximum inter-layer current 
density, which in turn can be used to probe vortex structure in the mixed state. 

We have carried out experiments on several k— (ET) 2 Cu(NCS)2 single crystals. The 
dimensions of each crystal were approximately 0.75x0.5x0.2 mm 3 , with the low conductivity 
axis the shortest of the three. We found that all samples gave qualitatively similar results. 
Microwave impedance measurements were carried out using a cavity perturbation technique 
described elsewhere We used the TEOln modes of cylindrical copper cavities (n = 

1, 2, 3, etc.), and the TE112 mode of a rectangular copper cavity. The use of a range of 
cavities enabled wide frequency coverage from 28 — 153 GHz. Typical loaded Quality factors 
for these cavities for a given measurement were approximately 5 x 10 3 to 2.5 x 10 4 . The 
samples were placed in one of two configurations within the cavity, which correspond to 



two different microwave field configurations |K|. The first position induces strictly in-plane 



currents, with the AC microwave magnetic field perpendicular to the layers (H^c -L be). 
The second configuration, with H^c || be, generates both in-plane and inter-layer Josephson 



currents; for a discussion of the electrodynamics of organic superconductors, see ref |TT . 
Both cavity positions were such that the DC magnetic field (B uc) was perpendicular to the 
conducting planes, i.e. B^c || a. Temperature control was achieved using a small resistive 
heater attached mechanically to the cavity, and a Cernox thermometer. For inter-layer 
measurements, DC magnetic fields were generated using an 8 T superconducting solenoid. 
The in-plane measurements were conducted in the 33 T resistive magnets at the National 
High Magnetic Field Laboratory in Tallahassee, FL. In both cases, field sweeps were made 
at a rate of approximately 1 T per minute. 

Figure 1 shows dissipation due to in-plane currents, plotted versus magnetic field, for 
temperatures in the range 1.6 to 4.5 K; the frequency is 44.4 GHz and the traces are offset for 
clarity. As expected for this configuration, the dissipation (oc surface resistance) increases 
monotonically with increasing DC field. Indeed, the surface resistance approximately follows 
a behavior, consistent with flux-flow type dissipation ||12|| . It is interesting to note the 
lack of (or extremely weak) temperature dependence in the various traces in Fig. 1 [13| . This 
is unexpected given that /i H c2 (~ 2 — 3 T at 2 K) varies by at least 0.5 T over the equivalent 
temperature range. At present, the origin of this behavior is not known and will form the 
basis for future investigations. The main point to note from the in-plane measurements is 
the lack of any resonant features (see Fig. 2 below). These findings are in stark contrast 
to recent studies by Schrama et al. [0], where it is claimed that an observed resonance is 



attributable to the in-plane response. 

The second configuration, in which a mixture of in-plane and inter-layer currents are 
exited {i.e. H^c || be), gives very different results from those shown in Fig. 1. Fig. 2 
plots dissipation versus magnetic field for several frequencies in the range 28 — 153 GHz; the 
temperature is 2.0 K in each case. A pronounced resonant feature can now be seen (indicated 
by arrows), which becomes sharper and moves to lower magnetic field upon increasing the 
frequency. Due to the absence of this resonance in the purely in-plane response (Fig. 1), 
we conclude that this resonance is related to dissipation within the sample caused by inter- 
layer currents [i.e. dissipation oc a a (u)}, and that insufficient care was taken to separate 



these contributions in the work of Schrama et al. JTj,|T5|. The anti-cyclotronic nature of the 
resonances shown in Fig. 2 has been well documented in the HTSs ||16|| , and has also been 
observed in previous work on this OS [|[, though only two frequencies were used. The wider 
frequency range employed here clearly confirms the anti-cyclotronic trend. 

Figures 3a and b show the temperature dependence of the resonance for uj p /2tt = 76 and 
111 GHz, respectively. Note that, as the temperature increases toward the superconducting 
to normal transition (T c ~ 10 K), the resonances vanish. This observation would seem to 
indicate that it is the superconducting carriers that are responsible for the sharp dissipative 
peak. This fact, along with the anti-cyclotronic nature of the resonance, and the AC current 
polarization dependence of the effect (i.e. resonances are not observed for Hac -L be), leads 
us to conclude that the resonance indeed corresponds to a Josephson plasma mode. Note, 
no apparent hysteresis was observed for the up and down sweeps of the magnetic field in 
this material, hence only up sweeps are plotted for clarity. 

Looking more closely at the 76 GHz data in Fig. 3a, we see that the resonance peak 
position decreases monotonically with respect to increasing temperature. This behavior is 
consistent with previous observations in this OS ||. However, we see an entirely different 
temperature dependence in the tail of the resonance — the structure in the low field tail (~ 
0.1 T) shows an initial increase in field with increasing temperature, reaching a maximum, 
then decreasing back to zero field as T c is approached. This so-called "cusp" behavior has 



been observed in some HTSs ||16|| , but never in any of the OSs. 

For the 111 GHz data in Fig. 3b, we see a similar resonance structure, but now the cusp 
behavior can be seen in the position of the resonance rather than its tail. The phenomenon 
responsible for the cusp has not changed, however, the JPR has been pushed to lower field by 
working at higher frequencies, thereby moving the resonance position into the cusp regime. 

Measurements spanning a broad frequency range, as shown in Fig. 4, indicate that the 
cusp behavior is confined to a relatively narrow field and temperature domain. In fact, 



comparison of the cusp position with published magnetization data JjTjTB show that it falls 



along, or very near to, the irreversibility line. Due to the change in nature of the temperature 



dependence of the JPR at the cusp position, we conclude that there is a transformation in 
the vortex structure at (or near) the irreversibility line. 

It is well known that a dimensional crossover phenomenon between a three-dimensional 
flux-line-lattice and a quasi-two-dimensional (Q2D) pancake vortex lattice occurs at approx- 
imately 10 niT in k— (ET) 2 Cu(NCS) 2 [ |T^ , [ZD| ] . Thus, it is assumed that, over a large portion 
of the superconducting phase diagram probed in this work (and shown in Fig. 4), there 
is strong disorder in the pancake positions as viewed along the a-axis (least conducting 
direction). At the same time, however, long range order may be present within the super- 
conducting layers (i.e. a Q2D vortex lattice) below some characteristic melting temperature 
[T m (S)]. We propose two possible explanations for the cusp phenomenon, both of which 
involve transitions in the vortex structure, and both of which account qualitatively for the 
observed B, T- dependence of the plasma frequency. 

One possible explanation for the cusp phenomenon involves a 2D solid to liquid melting 
transition pT| , p2| . This occurs when thermal (or quantum) fluctuations in the positions 
of the vortices become comparable to the inter- vortex separation, i.e. either when the 
flux density exceeds a critical value B m (T), or when the temperature exceeds a critical 
value T m (B). Below B m (T) and T m (B), in the Q2D pancake vortex lattice phase, random 
defects will tend to pin each Q2D solid in such a way that there is little correlation between 
the locations of vortices in adjacent layers, i.e. (cos <£> n) „ + i) will be small, and inter-layer 
Josephson tunneling will be suppressed. Upon increasing the temperature slightly, thermal 
fluctuations will tend to cause the vortices in each layer to deviate from their mean positions. 
Provided that these displacements are small compared to the distance between adjacent 
vortices, the effect of increasing temperature will actually result in an increase in the inter- 
layer tunneling probability and, hence, an increase in the maximum inter-layer Josephson 
current. Eventually, upon further increasing the temperature, the vortex lattice melts. In 
the liquid phase, increased thermal fluctuations (i.e. increasing temperature) will tend to 
suppress inter-layer tunneling and, hence, the maximum inter-layer Josephson current. Thus, 
completely opposite temperature dependences of the JPR frequencies are expected in the 
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Q2D solid and liquid phases (see Eq. 1), i.e. exactly as observed experimentally. 

Another possible explanation for the cusp phenomenon involves a depinning transition 
p3| , |2^] whereby, upon increasing the temperature, thermal fluctuations exceed some critical 
depinning threshold. As this threshold is approached, one can think of a "pinned" Q2D 
vortex lattice performing larger and larger collective displacements from equilibrium, until 
eventually it becomes completely "depinned," or mobile. In this sense, lattice melting and 
depinning represent qualitatively similar phenomena — indeed, under certain conditions, 



lattice melting and depinning may occur together [gj] . However, the physical origin of each 
transition is very different: in the case of melting, thermal fluctuations overcome the inter- 
vortex interactions which maintain long range Q2D order; in the case of depinning, thermal 
fluctuations suppress the influence of random defects, and long range Q2D ordering of the 
vortices persists. For both cases, one should expect opposing temperature dependences for 
(cos <fn,n+x) above and below the respective transitions, since it is pinning that is responsible 
for the suppression of interlayer tunneling at low temperatures, while the effects of this 
pinning are absent, or at least greatly diminished (see below), in the high temperature 
liquid or depinned states. 

We can now relate the observed magnetic field and temperature dependence of the res- 
onances to the above discussion. Application of a magnetic field tends to reduce the JPR 
frequency co> p (B,T) (see discussion below). Thus, for a fixed frequency measurement, a res- 
onance is observed when u; p (B,T) matches the measurement frequency. Consequently, any 
trend that increases u p (B,T) will shift the observed resonance to higher magnetic fields, 
since a stronger field is required to reduce u; p (B,T) so that it matches the measurement fre- 
quency. The converse is true for trends which decrease co> p (B,T). Thus, in the pinned Q2D 
solid phase, increasing temperature will tend to shift the resonance to higher field while, in 
the liquid or depinned Q2D solid phase, one expects the resonance to move to lower fields 
upon increasing the temperature. This is precisely the behavior observed in Fig. 4b, where 
the cusp marks the change in temperature dependence of the resonance. 

Next, we consider the frequency dependence of the resonance (at constant temperature). 



It has previously been assumed for this substance that the frequency of the resonance (u p ) 
obeys a power-law field dependence of the form 



v 



A5^exp(T/T c ), (2) 



where A is a constant. The inset to Fig. 5 shows two fits to a power-law for T = 2.0 K: the 
dashed line is a fit to data covering a wide field range; the dotted line is a fit to data above the 
melting/depinning transition only. The former fit gives a value for the exponent of fi = 0.7, 
which is consistent with previous work where only two data points were available [|]]. The 
latter fit gives a value for the exponent of /i = 1.7. Overall, the power-law does not appear 
to be a good fit, even when restricted to data above the melting/depinning transition, i.e. 
closer inspection of the inset to Fig. 5 shows that the exponent \x increases with increasing 
field. We also note that an exponent /i > 1 is unphysical based on the various theories 
which predict a power-law field dependence of the squared plasma frequency P,p5|; \i of 
order unity corresponds to little or no correlation between vortices in adjacent layers, while 
H < 1 indicates some degree of ordering of the vortices along lines parallel to the a-axis. In 
the limiting case of a vortex liquid, fj, — 1 has been predicted f2"5| , and /i close to unity has 
been observed experimentally in BSCCO 

In the main panel of Fig. 5, we show two further fits to the data using a week pinning 
theory developed by Bulaevskii et al. ||, where a decaying exponential field dependence of 
ajp(B,T) is predicted; the dashed and dotted curves correspond to the same field ranges as 
the power-law fits discussed above. In this case, 

ul(BJ) = ^J exp[-(B/B D )\ 

where (3) 
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Bd is a de-coupling field, above which inter-layer phase coherent tunneling ceases f26|l ; 7 Q 



is the anisotropy parameter; E p is a pinning parameter; and An is the in-plane penetration 
depth [|9] . Clearly, the above expression yields superior fits to the data over both field ranges 
(there is only a ~ 10% difference in the fitting parameters obtained from the dashed and 
dotted curves). 

At the present time it is not clear why the Bulaevskii theory, which was developed for 
a weakly pinned vortex solid, fits our data so well. In particular, since most of the data 
in Fig. 5 were obtained above the melting/depinning transition, one might expect a 1/B 
dependence of co> 2 , as predicted for a vortex liquid phase [25fl. One possible conclusion is 



that quasi-long range 2D order persists in this high temperature regime, and that pinning 
continues to play a role. Further detailed frequency dependent studies will be required to 
resolve this issue. 

It is interesting to note that, if one adopts the weak pinning theory, it is possible to 
extract certain sample parameters for comparison with published results. Using e c = 25, 7 Q 
= 100—200, and An =1 /im [0], we find the zero field maximum interlayer current density, 
J Q = 2 x 10 3 A cm -2 . From this, we are able to deduce the interlayer penetration depth, 
X± = 90 //m, which is on the order of previous measurements [0]. Further, we find a value 
for the decoupling field, = 0.7 T, from which we can extract the pinning parameter E p 
= 6.3 x 10 -13 to 1.3 x 10~ 12 J 2 m -4 . This range of values is three orders of magnitude lower 



than the BSCCO measurements discussed in ref. ||27|| . This is not surprising, since it is well 
established that the ET salts are considerably cleaner systems than the HTSs, thus having 
far fewer pinning sites. Moreover, for the weak pinning limit, it has been noted that the 



in-plane critical current depends on the pinning parameter given above [27 



J o6 = {±c/m o ){V p /?>) x l\ (4) 

Using the range of values above for E p , we find the critical in-plane current density to be in 
the range, J a b = 3 x 10 4 A cm~ 2 to 4 x 10 4 A cm -2 , consistent with values obtain by other 



methods |4j]. Thus our measurements have proven consistent with previous work, and seem 
to agree excellently with the weak pinning theory of Bulaevskii et al. || . 

Finally, we have noted a slight sample dependence in our measurements, as well as some 
variance compared to other published results PJI^|. The most notable differences can be 
found in the resonance field positions for a given frequency and temperature, as is clearly 
seen in Fig. 4, where several of the curves corresponding to different frequencies lie on top 
of each other and even cross in one case. This sample dependence is, perhaps, not surprising 
given the dependence of u; p (B,T) on the pinning parameter, E p , which is liable to depend 
on sample quality. However, we also see an apparent scatter in the cusp positions in Fig. 3. 
At this stage, it is unclear whether this scatter is real, i.e. attributable to sample quality, 
or whether it is associated with uncertainties introduced during data analysis. Establishing 
whether there is a sample dependence in the cusp positions could help resolve the nature 
of the proposed transition in the vortex structure/dynamics; a first order melting transition 
should not be influenced by sample quality, whereas one would expect a depinning transition 
to exhibit pronounced sample dependence (see refs. p|,pE|, for recent results on the K-phase 
ET salts). 

In conclusion, Josephson plasma resonance measurements have provided a wealth of 
information regarding the superconducting properties of the organic charge transfer salt 
k— (ET) 2 Cu(NCS) 2 . Due to the extremely clean nature of this material, it is apparent that 
the power-law field dependence of the squared plasma frequency, as observed in several 
High— T c compounds, is inapplicable in the present case. Instead, we are able to fit our data 
with a weak pinning theory, from which we obtain excellent agreement with several well 
known material parameters, e.g. the inter-layer penetration depth, and the in-plane critical 
current. This novel technique also allows us to detect subtle changes in the Q2D vortex 
structure/dynamics. In particular, we see a cusp in the temperature dependence of the 
plasma frequency u; p (B,T) which may be due to either vortex lattice melting or depinning. 
Further examination of this cusp should lead to a better understanding of the observed 
transition, and to the properties of vortex matter in general. 
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Figure captions 

Fig. 1. Field dependent dissipation (oc surface resistance) due to in-plane currents, at 
various temperatures, and for uo/2-n = AAA GHz (offset for clarity). As expected for flux 
flow type dissipation, the surface resistance varies approximately as B 1 / 2 for B< /i H c2 . This 
configuration shows no resonant behavior and little temperature dependence. The inset 
depicts the geometry of the crystal, the orientations of the applied DC and AC magnetic 
fields, and the AC current paths around the edges of the sample. 

Fig. 2. Magnetic field sweeps, at various frequencies, with H^c || be; the temperature 
is 2.0 K in each case. Notice the sharp resonance feature, indicated by arrows, which is 
attributable to inter-layer Josephson currents. Note that the resonance field position, B G , 
decreases with an increasing frequency, i. e. anti-cyclotronic behavior. 

Fig. 3. Temperature dependence of the Josephson plasma mode for a) uj/2ti = 76 GHz and 
b) uo/2n = 111 GHz. Notice that in a), the resonance peak position decreases monotonically 
as a function of temperature; however, the structure in the tail of this resonance exhibits a 
"cusp" behavior. In contrast, due to the anti-cyclotronic behavior of the JPR, the resonance 
in b) has moved to lower field, and we now see a cusp behavior in the resonance position. 
The inset depicts the geometry of the crystal, the orientations of the applied DC and AC 
magnetic fields, and the AC current paths across the surfaces of the sample. 

Fig. 4. a) Plot of the cusp positions in field, versus temperature, for several samples and 
frequencies; the inset plots several representative series of resonance fields (for a single sam- 
ple and a range of frequencies) versus temperature, and the method used to determine the 
location of the cusp, i.e. the crossing points of the extrapolated fits (dashed lines) to the 
low and high temperature behavior. The cusp positions in the main panel of a) are scat- 
tered around the irreversibility line (solid line), b) Resonance field B c vs. temperature for 
measurements spanning a wide frequency range (28-153 GHz), and for several samples. The 
higher frequencies exhibit a cusp behavior, while lower frequencies show a monotonic tem- 
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perature dependence. The dashed line indicates the approximate position of irreversibility 
line, and the solid lines are guides to the eye. 

Fig. 5. When <Ji is plotted against B 3//2 , an exponential decay is observed. This behavior 
is predicted for a weakly pinned vortex lattice || . The dashed line is a fit to all of the data 
points, while the dotted line represents a fit to data above the melting/ depinning transition. 
Inset: we have tried to fit the data with a power-law, as has been observed in the HTSs, and 
previously claimed in this material [3]. However, this is clearly not a good fit when a wide 
field/frequency interval is sampled. The dashed line represents a fit to all of the data, giving 
an exponent \i = 0.7, while the dotted line is a fit to data above the melting/depinning 
transition, giving /i = 1.7. 
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FIG. 1. M. Mola et al. 
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FIG. 2. M. Mola et al. 
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FIG. 3. M. Mola et al. 
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FIG. 4. M. Mola et al. 
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FIG. 5. M. Mola et al. 
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